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The formation of diatomic molecules with rotational and vibrational coherence is demonstrated 
experimentally in free-to-bound two-photon femtosecond photoassociation of hot atoms. In a ther- 
mal gas at a temperature of 1000 K, pairs of magnesium atoms, colliding in their electronic ground 
state, are excited into coherent superpositions of bound rovibrational levels in an electronically 
excited state. The rovibrational coherence is probed by a time-delayed third photon, resulting in 
quantum beats in the UV fluorescence. A comprehensive theoretical model based on ab initio cal- 
culations rationalizes the generation of coherence by Franck-Condon filtering of collision energies 
and partial waves, quantifying it in terms of an increase in quantum purity of the thermal ensemble. 
Our results open the way to coherent control of a binary reaction. 
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When a pair of colliding atoms is irradiated by a fem- 
tosecond pulse, an important scenario is photoassocia- 
tion (PA) in which a chemical bond is "photo-catalysed" 
via a free-to-bound broadband optical transition [lHa]- 
Femtosecond laser pulses can be shaped in their ampli- 
tude, phase and polarization. They are thus the basic 
tool for coherent control which employs interference of 
matter waves to steer the dynamics of a quantum system 
toward a desired outcome 0-Q- Chemical reaction con- 
trol, as a means for photo-inducing new chemical synthe- 
sis, is a long-standing yet unrealized goal from the early 
days of coherent control [8j . While coherent control has 
been extremely successful for destructive unimolecular 
processes such as ionization or dissociation 10l4l4| , con- 
trol of associative binary reactions still remains an un- 
resolved puzzle. A controlled photoassociation step be- 
tween two atoms or two molecules has been proposed as 
a first step in an extended control scheme of more com- 
plicated bimolecular reactions [15J. Femtosecond pho- 
toassociation leading to chemical bond formation is thus 
the prerequisite for coherent control of a binary reaction. 
Femtosecond photoassociation of hot atoms (hot fs-PA) 
has been demonstrated in a singlepioneering experimen- 
tal study more than 15 years ago |l[, using mercury and a 
single-photon UV transition. Coherent rotational molec- 
ular dynamics were observed following the PA. In order to 
utilize hot fs-PA as a basis for coherent control of binary 
reactions, it is, however, essential to induce vibrational 
coherence of the photoassociated molecules since the vi- 
brational degrees of freedom determine the fate of bond 
making and breaking. In view of a large control toolbox, 
it would be optimal to have both one-photon and multi- 
photon fs-PA available. The high peak powers of fem- 
tosecond laser pulses can easily drive multi-photon tran- 
sitions, resulting in a large range of molecular excitation 



energies amenable to fs-PA. Multi-photon PA allows for 
accessing electronically excited molecular states that are 
inaccessible in one-photon PA due to the different selec- 
tion rules. It also allows for utilizing multi-photon control 
strategies that differ from one-photon control strategies 
due to an extended manifold of state-to-stat e q uantum 
pathways and large non-resonant Stark shifts [16| . More- 
over, in view of experimental feasibility, fs-PA using non- 
resonant multi-photon transitions with near-IR or visible 
photons is desirable in order to associate molecules which 
have their lowest electronic transition in the deep-UV or 
VUV spectral range where pulse shaping is more difficult. 
Multi-photon fs-PA thus extends the variety of molecu- 
lar species and reactions that are candidates for coherent 
control using shaped pulses. 

In this Letter, we demonstrate two-photon hot fs- 
PA with subsequent coherent molecular motion that is 
of both rotational and vibrational nature. Our exper- 
imental and theoretical results are presented for fs-PA 
of hot magnesium atoms into bound magnesium dimer 
molecules, i.e., Mg+Mg— ^Mgj, with the thermal ensem- 
ble of Mg atoms held at a temperature of 1000 K. We 
show that femtosecond laser pulses can carve from an 
initial incoherent mixture a sub-ensemble with increased 
quantum purity, Tr[p ], and dynamical coherences that 
are amenable to coherent control. 

Magnesium is experimentally well-suited for PA stud- 
ies for two main reasons, (i) The experimental back- 
ground signal originating from direct bound-bound op- 
tical excitation of already-bound ground-state dimers is 
very small. This is due to the van der Waals nature 
of the molecular ground electronic state whose shallow 
well of only about 400 cm -1 is outweighted by the ro- 
tational barrier and does not support any bound states 
for most of the thermally populated partial waves, (ii) 
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FIG. 1: The excitation scheme of the two-photon femtoseo 
photoassociation of magnesium. The shaded areas indicate 
the thermally populated scattering states in the electronic 
ground state and the vibrational band of IT 9 -state molecules 
inferred from the autocorrelation function. The position of 
the vibrational band above the weak-field two-photon reso- 
nance is rationalized by the strong dynamical Stark shift. 
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FIG. 2: (color online) Two-photon femtosecond photoassocia- 
tion of hot magnesium atoms: (a) Experimental pump-probe 
signal showing the measured UV fluorescence intensity as a 
function of the pump-probe delay (black line) . The red line 
corresponds to rotational dynamics only. It results from a low- 
pass step filter set to include only frequencies smaller than 
4THz or 133.42 cm -1 . (b) Fourier spectrum of the experi- 
mental pump-probe trace (solid black and grey lines) . Super- 
imposed are the theoretical amplitude spectra corresponding 
to purely vibrational dynamics (blue dotted line), purely ro- 
tational dynamics (dashed red line) and vibrational dynamics 
including rotational broadening (purple dotted line). 



The electronically excited states of Mg2 with deep poten- 
tial wells are conveniently accessible in a free-to-bound 
PA process via optical transitions of two or three NIR 
photons. Figure Q] displays the relevant Mg 2 potential 
energy curves and the present excitation scheme. The 
inset of the figure shows the ground -X^E* state for 
three different values of partial waves, J, of the col- 
liding Mg atoms: Already for J w 50, the X-state 
does not support any bound levels. The temperature 
of T — 1000 K corresponds to a collision energy UbT/2 
of ~ 350 cm -1 . Pairs of Mg atoms colliding in the A" 1 E+ 
electronic ground state of Mg 2 are excited via a broad- 
band free-to-bound non-resonant two-photon transition 
into coherent superpositions of bound rovibrational lev- 
els in the electronically excited (l) 1 !^ state. The non- 
resonant two-photon coupling is provided by off-resonant 
dipole transitions to intermediate states of ungerade sym- 
metry with the largest contribution coming from the low- 
est states, A X E+, (l) 1 !^ and ^^II,,. PA is induced by 
pulses of 70 fs transform-limited duration, 840 nm central 
wavelength, linear polarization, and a transform-limited 
peak intensity of about 5xl0 12 W/cm 2 . A time-delayed 
probe pulse, identical to the pump pulse, probes the dy- 
namics of the excited dimers via a one-photon excitation 
to the higher lying (l) 1 !^ electronic state. This state 
has a strong one-photon allowed transition to the A X E+ 
ground state with a lifetime of a few ns. The intensity of 
the resulting UV fluorescence at 285-292 nm is measured 
as a function of the pump-probe time delay r, yielding 
the pump-probe signal transient. The experiments are 
conducted in a static cell with argon buffer gas and a 



pressure of the Mg vapor of ~ 5 Torr. The sample is ir- 
radiated at a repetition rate of 1 kHz by the pump and 
probe laser beams. UV fluorescence emitted toward the 
laser-beam entrance to the cell is collected at a small an- 
gle from the laser-beam axis using an appropriate optical 
setup. 

Figure [5] presents experimental results for two-photon 
fs-PA of hot magnesium atoms, showing the measured 
pump-probe transient in Fig. [2ja) and its Fourier Trans- 
form (FT) spectrum in Fig. &b). The time-dependent 
pump-probe transient signal, symmetric around r = 
since pump and probe pulses are identical, is ob- 
served also for delay times where the photoassociating 
pump pulse has already ended. As corroborated by the 
calculations presented below, observation of this time- 
dependence demonstrates that (i) molecules have been 
formed in a non-stationary coherent state, i.e., in a coher- 
ent superposition of molecular rovibrational states, and 
(ii) the generated ensemble of photoassociated molecules 
must be of increased purity as compared to the initial 
thermal, incoherent atomic ensemble - without an in- 
crease in purity, no coherence could be observed. The 
FT spectrum of the experimental pump-probe signal, cf. 
Fig. (Hb), is composed of a narrow low- frequency peak 
at around 38.9 cm' 1 and a broad high-frequency peak 
located at 189.6 cm -1 . It shows a zero spectral phase, 
confirming that pump and probe excitation occur at the 
same R, cf. Fig. Q] Below, we show the low-frequency 
peak to be associated to coherent rotational dynamics 
and the high-frequency peak to coherent vibrational dy- 
namics. 

To interpret the experimental data, we have con- 



structed a comprehensive theoretical model. State-of- 
the-art ab initio methods have been employed to calcu- 
late the Mg 2 potential energy curves presented in Fig. Q] 
as well as the corresponding one-photon and two-photon 
transition dipole moments, dynamical Stark shifts, and 
spin-orbit interaction, non-adiabatic and radial couplings 
in the electronically excited states. The short-range part 
of the potentials and transition moments were calculated 
within the framework of the equation of motion (reponse) 
coupled cluster method restricted to single and double 
excitations [17|, |l8| in a large basis set. The long-range 
part of the potentials as well as the non-adiabatic, ra- 
dial and spin-orbit couplings have been obtained using 
multi-reference configuration interaction. A high accu- 
racy of the computed potential energy curves is con- 
firmed by a good agreement of the calculated spectro- 
scopic constants with their experimental values. Our ab 
initio data are by far much more accurate than previ- 
ous results reported in Refs. [19(. The PA dynamics 
are simulated non-perturbatively by solving the time- 
dependent Schrodinger equation on a multi-surface grid 
including all non-adiabatic couplings and one-photon and 
two-photon transition dipole operators. It is crucial to 
correctly represent the initial thermal ensemble of scat- 
tering atoms (as opposed to a single state which is al- 
ways coherent). The incoherent initial state is described 
by the density operator in the canonical ensemble, p(t — 
0) = cxp[-/3H 3 ]/Z with /3 = l/k B T, H g the Hamiltonian 
for nuclear motion in the ground electronic state, and 
Z the partition function, Z — Tr[exp[— /3H g }] It is con- 
structed using random-phase wavefunctions [20]. Details 
of the electronic structure calculations and the random- 
phase thermal wavefunction approach are reported else- 
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Separating rotational and vibrational motion, the PA 
dynamics is simulated by propagating many realizations 
of random-phase thermal wavefunctions under the pump 
and probe pulses for each value, J, of the incoming par- 
tial waves. Computing the expectation value of an ob- 
servable for each propagated random-phase wavefunc- 
tion and incoherently summing all individual expecta- 
tion values yields the thermally averaged time-dependent 
expectation value. For example, the calculated spec- 
trum shown in Fig. [U[b) (blue dotted line) is obtained 
by Fourier transforming the thermally averaged time- 
dependent population of the 1 IL U states. Since the full 
pump-probe simulation for many random-phase realiza- 
tions and all thermally populated partial waves J is 
numerically rather demanding, our calculations account 
only for optical transitions of AJ = 0. Our quantum 
dynamical results therefore reflect the dependence of hot 
fs-PA on the partial waves, J, but the calculated dynam- 
ics is of purely vibrational nature. 

Comparing the experimental spectrum to the calcu- 
lated one (black solid and blue dashed lines in Fig.[2](b)), 
we find the high-frequency peak to be due to coherent 



vibrational dynamics. The range of populated vibra- 
tional levels in the ^g state is indicated by the grey 
shaded area in Fig. [1] We approximately account for 
the rotational dynamics by carrying out a quantum- 
beat analysis: We assume rovibrational 1 II 5 levels with 
J' = J, J ± 1, J ± 2, J ± 3, J ± 4 to get excited by the 
pump pulse from each incoming partial wave J. Each 
pair of rovibrational 1 FI ff levels, (v[,J[), (v' 2 ,J 2 ), with 
v' 2 = v' 1: v[ + 1 and \J 2 — J[\ = 2,4 contributes to the 
FT spectrum an amplitude Pin (J) at the frequency 
wqb = (Ev'^Ji ~ K' 2 ,J0/K- p in g (J) is the 1 II a state 
population obtained in the vibrational dynamics simula- 
tion for each incoming J, shown in Fig. [3(a). The re- 
sult of the quantum-beat analysis is shown by the red 
dashed line in Fig. EJb) and describes the low-frequency 
peak in the experimental spectrum very well. We thus 
find this peak to fit pairs of 1 H g levels with v[ = v 2 
and J{, J 2 ~ 40 . . . 100. In particular, the quantum beat 
frequency corresponding to the maximum of the peak 
matches only pairs with J{ 2 ~ 50 . . . 60. Combining our 
approximate treatment of the rotational dynamics with 
the full quantum simulation of the vibrational dynamics 
leads to rotational broadening of the vibrational peak, cf. 
purple dotted line in Fig. [2jb) . We attribute the remain- 
ing small discrepancy between the experimental and the 
theoretical spectrum to the separation of rotational and 
vibrational dynamics and the uncertainty of the ab initio 
calculations. In particular, the sharp feature in the cal- 
culated vibrational spectrum is due to the non-adiabatic 
couplings between the ^h states which carry the largest 
error bar of the otherwise very accurate ab initio results. 
To summarize, the peaks observed in the experimental 
spectrum originate from coherent excitation and probing 
of pairs of rovibrational levels (v[, J{) and (i> 2 , J 2 ) with 
high rotational quantum numbers in the (l) 1 II g electron- 
ically excited state. 

As illustrated in the inset of Fig. [TJ for values of J > 50, 
the X 1 ^ ground state of Mg 2 does not support any 
bound molecular level. We therefore deduce from our 
quantum beat analysis that indeed the majority of the 
bound excited Mg 2 molecules, that we probe experimen- 
tally, have been generated in a free-to-bound coherent PA 
process. This conclusion of free-to-bound PA of high- J 
molecules is further supported by the calculated relative 
PA probabilities, shown in Figure [3]( a), as a function of 
the incoming partial wave, J. Maximum PA probability 
is observed for J = 55, for which the electronic ground 
state does not support bound levels, cf. Fig.[TJ A careful 
analysis shows that about 80 per cent of the total signal 



is due to free-to-bound PA [21 1 



In order to rationalize and quantify the generation 
of the coherence, that is observed in Fig. [2] out of 
a completely incoherent initial ensemble, we consider 
the change in quantum purity, Tr[/3 ], due to the fem- 
tosecond laser pulses. The initial purity, determined 
by the temperature, T — 1000 K, and density, p = 
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FIG. 3: (a) Calculated relative PA probability as a function of 
incoming partial wave J. (b) Quantum purity, Tr[p(tf ina i) 2 ] 
and dynamical coherence of the sub-ensemble of molecules in 
the (1) 1 II 9 state vs laser intensity 



4.8 • 10 16 atoms/cm 3 , of the experiment, is estimated in 
terms of the purity of a single atom pair represented in 
our computation and the probability, P2, of finding two 
atoms in our computation volume, V. Given p and V = 
4/37rr 3 = 6.33-10~ 21 cm 3 , we find p 2 = 9.2-1CT 8 . The pu- 
rity of a single atom pair in our computation box is given 
by the vibrational purity weighted by the probability, 
Pj , for occupation of the sub-ensemble with angular mo- 
mentum J. For an atom pair in their electronic ground 
state, the purity becomes V 9 — "p\ ^2j PjVj where 

V 9 j = Tr[(e-^) 2 ]/Zj, and H^ and Zj are the Hamil- 
tonian for vibrational motion and partition function, re- 
spectively, of partial wave J. Evaluating Vj with thermal 
random-phase wavefunctions, we obtain V 9 = 2.6 • 10~ 20 
for the initial purity. For the ensemble of molecules in 
the electronically excited 1 H g state, the density opera- 
tor is given by p e — J2jPjP"j- Here, Pj is the prob- 
ability for occupation of the excited state sub-ensemble 



with angular momentum, J, cf. Fig. EJa). pj is the 
normalized density operator of the excited state J sub- 
ensemble which is constructed by incoherently averag- 
ing dyadic products of the propagated thermal wavefunc- 
tions, \i>n,j(tf ina i))(i>n,j(t final)], over all random phases 
n and partial waves J {t final denotes the time when the 
pump pulse is over) . The purity for the excited state sub- 
ensemble is then given by V e = J2j(Pj) 2 Tr[(/5}) 2 ]. We 
obtain a purity P e — 3.04 • 10~ 4 for the molecular sub- 
ensemble in the 1 II ff excited state for the experimental 
pulse parameters. We thus observe a dramatic increase in 
the quantum purity, Tr[/3 ], induced by the femtosecond 
laser pulse. The underlying physical mechanism can be 
viewed as " Franck-Condon filtering" : for a given initial 
J value there is only a limited range of collision energies 
that allow the colliding pair to reach the Franck-Condon 
window for PA located at short internuclear distances 121 ■ 



In order to obtain a quantitative estimate of the de- 
gree of distillation achieved by the fs-PA process, we 
have calculated the purity of ensemble of photoassoci- 
ated molecules in the 1 FI ff state for a range of laser in- 
tensities, cf. Fig. Elb). For weak fields, the purity is 
roughly constant as a function of intensity and almost 
three times larger than the purity obtained for the in- 
tensity of 5 x 10 12 W/cm 2 used in the experiment. As 
intensity is increased, a sudden drop in the purity is 
observed which levels off at large intensities. We at- 
tribute this drop to power broadening for strong fields, 
which brings more atom pairs into the Franck-Condon 
window for PA. The purity increase with respect to the 
initial ensemble is nevertheless many orders of magnitude 
for both weak and strong field. To further analyze the 
generation of quantum coherence, it is possible to sepa- 
rate static and dynamic contributions, p = p stat + Pd yn - 
This is most easily achieved in the energy representation 
where the static (dynamic) part corresponds to the di- 
agonal (off-diagonal) matrix elements. The dynamical 
contributions are quantified by the coherence measure 
C = Tr[p^ yn (t)] < Tr[p 2 (i)] Q. Tracing the coherence 
measure of the excited state, C e , as a function of laser 
intensity (red open diamonds in Fig. E)(b)) shows that 
most of the purity comes from the dynamical contribu- 
tion (C e = 2.86 • 1CT 4 for / = 5 x 10 12 W/cm 2 ). 

In summary, we have demonstrated femtosecond pho- 
toassociation of hot atoms using non-resonant two- 
photon transitions driven by NIR light. Quantum beats 
in the observed UV fluorescence reflect both vibra- 
tional and rotational coherence of the photoassociated 
molecules. Our experiment is the first to demonstrate 
generation of vibrational coherence in a light-induced bi- 
nary reaction. This is rationalized by our theoretical 
model in terms of Franck-Condon filtering, yielding an 
increased quantum purity of the sub-ensemble of the re- 
action products, the photoassociated molecules, as com- 
pared to that of the reaction partners, the atoms. Vi- 
brational coherence is a prerequisite for further control 
of the reaction products. Our work thus opens the way 
toward coherent control of a binary reaction. 
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